一个解析三维非结构网格生成方法

给出了一个三维非结构网格生成方法。提出弱区域指示函数描述三维封闭流场形状,并证明了该函数的存在不唯一性。利用该函数基于Delaunay方法依据流场尺度生成网格结构。其中,边界引力有效解决了流场边界恢复和内嵌边界问题,弹簧振子和正四面体趋进技术有效调整网格结构并提高网格质量。最后,弹簧振子和边界引力随着弱区域指示函数的演变与监测函数的指导成功应用到自适应网格生成。     

共轴刚性旋翼悬停状态气动干扰机理

建立了一个基于Navier-Stokes(N-S)方程的共轴刚性旋翼气动干扰数值模拟方法。应用运动嵌套网格技术模拟双旋翼反转运动。通过与试验值对比,验证了方法的有效性。分析了共轴刚性旋翼悬停状态的气动性能和流场特征,结果表明,双旋翼气动干扰主要来自4个方面:双旋翼尾迹涡相互诱导引起"涡诱导效应",使上旋翼气动性能优于下旋翼;双旋翼周期性相遇-离开过程中桨叶附着涡干扰引起"载荷效应",对应拉力周期性升降波动;双旋翼相遇时"厚度效应"使双旋翼拉力产生相反的脉冲波动;上旋翼尾迹涡与下旋翼桨叶碰撞引起垂直"桨-涡干扰效应",使下旋翼桨叶展向拉力分布受到干扰。     

楔形体入波浪水面数值模拟

数值模拟分析楔形体入波浪水面流体力学现象和运动姿态变化历程，为水上飞机在波浪水面起降提供技术和理论支持。基于有限体积法求解非定常雷诺平均Navier-Stokes（URANS）方程和SST k-ω湍流模型，结合整体动网格法的速度入口造波方法、阻尼消波法和流体体积（VOF）液面捕捉技术构造数值波浪水池，耦合三/六自由度模型，实现楔形体入线性规则波和不规则波的数值模拟。在楔形体入静水面算例验证的基础上，详细研究二维楔形体在线性规则波波峰、波谷、平衡位置-上行速度和下行速度处入水过程的受力特性和运动姿态变化。结果表明：速度入口造波方法的波浪数值解与理论解析解较吻合，偏差为1%；在二维楔形体入规则波过程中，楔形体所受垂向力和垂向速度变化趋势相同，对横向位移影响较小，其数值变化少于0.01；在平衡位置入水过程中，模型的滚转角和横向速度变化明显，其数值变化约为波峰、波谷位置处的8倍和10倍。分析原因：平衡位置处波浪内流速度变化较快，模型所受波浪力冲量较大，且楔形体两侧斜边与波浪的相互作用力也参与到模型的运动中。数值模拟二维楔形体在线性不规则波5个随机时间点入水及三维楔形体在波峰、平衡位置-上行速度处入水过程，其变化规律及形成原因与入规则波相似。     

复杂外形网格生成技术

复杂外形网格生成是进行流场数值模拟的一个关键技术及前提。本文论述了采用保角变换方法，求解偏微分方程方法以及它们在代数方法结合生成二维及三维具有多块网格结构的贴体且与边界正交网格的技术。     

用任意空间曲线控制的复杂组合曲面数控加工算法研究

提出了一种基于参数曲面三角化的复杂曲面刀具轨迹生成方法。即利用给定的任意空间曲线作为控制边界，构造新的参数曲面，在此简单参数曲面上生成刀具轨迹。离散化复杂组合曲面为若干平面三角片，投影已生成的刀具轨迹到平面三角片上，再考虑加工精度、余量、干涉等，便可得到复杂组合曲面的刀具轨迹。根据该方法还可以获得指定加工区域的刀轨。     

高速平板着水数值模拟

探索和揭示物体入水冲击的流体力学现象与机理对飞行器水上迫降问题的研究有重要的参考价值。对高速平板着水涉及到的复杂物理问题展开数值模拟,采用有限体积法求解非定常雷诺平均Navier-Stokes(URANS)方程和标准k-ε湍流模型,流体体积(VOF)模型捕捉水气交界面,整体动网格技术处理平板与水面的相对运动。在二维楔形体入水冲击的算例验证基础上,详细研究平板高速着水引起流体喷溅、射流、空气垫等现象和平板底面压力变化历程,结果表明:空气垫现象明显,俯仰角4°平板下表面出现规律的空气泡,10°时则不存在;平板下表面的水体沿壁面运动,当俯仰角为10°时,壁面水体的运动速度显著增加;在大俯仰角的情况下明显出现负压区。     

关键设计参数对摆线桨气动性能影响

以西北工业大学自行研制的摆线桨飞行器为研究对象,对简化的二维摆线桨模型进行了非定常数值模拟.在数值模拟模块中,桨叶的公转及俯仰振荡运动采用弹簧近似光滑模型和局部重划模型相结合的动网格技术来处理.重点研究了关键设计参数对摆线桨气动性能的影响,结果表明:随着桨叶数的增加,悬停气动效率提高;随着翼型厚度增加,推力变大,气动功耗减小;桨叶俯仰轴位置位于桨叶弦向中部位置时的功率载荷最大,悬停气动效率最高;随着最大俯仰角增大,气动功耗逐渐增加,悬停气动效率降低;当桨叶上、下半周最大俯仰角之和一定时,采用上半周最大俯仰角小的设置时,推力和气动功耗较大,悬停气动效率也更高.      

基于混合笛卡儿网格方法的可压流动数值模拟

以可压缩黏性流动的数值模拟为研究背景,发展了一套自适应混合笛卡儿网格(AHCG)方法以及基于有限体积方法的雷诺数平均Navier-Stokes(RANS)的数值求解方法.为更好地模拟边界层的黏性流动在近壁面处采用贴体结构网格,剩余计算区域自动生成与之相重叠的笛卡儿网格,并同时发展了基于流场特征的笛卡儿网格自适应技术.结合ADT(alternating digital tree)算法显著减少了网格生成中“挖洞”和“贡献单元”搜索的消耗机时,50万左右的网格数目下,搜索耗时为0.062s,仅为普通遍历方法的1/1847.通过二维圆柱与两段翼型绕流的数值算例显示,定常AHCG方法能够准确地预测物面压力分布与升阻力系数并且具备处理复杂外形的能力;同时通过二维非定常圆柱绕流问题与NACA0015矩形机翼翼尖尾涡的捕捉算例显示,结合了动态自适应网格加密的非定常AHCG方法尤其适用于旋涡主导流动.      

主要参数对面齿轮传动噪声辐射的影响分析

为研究面齿轮传动力学特性及噪声辐射,通过建立面齿轮传动的弯扭耦合动力学模型,在分析轴承支承处动载荷的基础上,结合有限元法和边界元法,形成了面齿轮传动的噪声辐射计算流程.此外,利用该计算流程,开展了算例分析,讨论了单一几何参数和工况对面齿轮传动噪声辐射的影响规律.仿真结果表明:模数和传动比对改善面齿轮传动噪声辐射有益,压力角对其噪声辐射影响不明显.研究结果为低噪声面齿轮传动的设计奠定了必要的理论基础.      

基于刚性动网格技术的动导数数值模拟

基于刚性动网格技术,建立了俯仰动导数的非定常数值计算方法.首先使用小幅度强迫俯仰振荡方法求解俯仰组合动导数,然后利用小幅度强迫升沉振荡方法求解洗流时差导数,通过两者相减即可得到俯仰阻尼导数.利用国际动导数标准模型Finner导弹进行算例验证,计算得到的俯仰组合动导数与试验值误差为2.76%,洗流时差导数值约为俯仰阻尼导数的11.5%,与文献的结果一致.结论表明:动导数单独模拟方法具有较好的工程实用价值,且可以推广到横向以及航向的动导数数值模拟.